Finite-Size Effects in Surface-Enhanced Raman Scattering from Molecules Adsorbed 

on Noble-Metal Nanoparticles 
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We study the role of strong electron confinement in surface-enhanced Raman scattering from 
molecules adsorbed on small noble-metal particles. We describe a new source of Raman signal 
enhancement which originates from different behavior of sp-band and d-band electron densities 
near the particle boundary. In small particles, a spillover of sp-electron wave- functions beyond the 
classical radius gives rise to a thin layer with diminished population of d-electrons. In this surface 
layer, the screening of sp-electrons by d-band electron background is reduced. We demonstrate 
that the interplay between finite-size and underscreening effects results in an increase of the surface 
plasmon local field acting on a molecule located in a close proximity to the particle boundary. Our 
calculations, based on two-region model, show that the additional enhancement of Raman signal 
gets stronger for smaller nanoparticles due to a larger volume fraction of underscreened region. 

PACS numbers: 33.20.Fb, 78.67.Bf, 71.45.Gm, 33.50.-j 



I. INTRODUCTION 

A renewed interest in surface-enhanced Raman scat- 
tering (SERS)P, 01 stems from the discovery of an enor- 
mous (up to 10^^) enhancement of single-molecule Ra- 
man signal in silver nanoparticle aggregates 

ii Al- 

though the relative importance of various mechanisms of 
SERS is still an issue under active discussion, the main 
source is attributed to the electromagnetic enhancement 
due to the local field of surface plasmon (SP) excited 
in a nanoparticle by the incident light 0, |a S 01 (see 
Fig. nj. Other possible enhancement mechanisms in- 
volve dynamical charge transfer between a nanoparticle 
and a molecule (chemical mechanism) and have been ad- 
dressed , e.g., i n Refs. 0, 0, 1^ ■ Recent experimen- 
tal [llfll [l5lllllT7t and theoretical [TlligllM I2ll23 | 
studies indicate that the anomalously strong Raman sig- 
nal originates from "hot spots" - spatial regions where 
clusters of several closely-spaced nanoparticles are con- 
centrated in a small volume. The high-intensity SERS 
then originates from the mutual enhancement of SP local 
electric fields of several nanoparticles that determine the 
dipole moment of a molecule trapped in a gap between 
metal surfaces. 

Although in single nanoparticles the magnitude of 
SERS is considerably smaller, with enhancement up to 
10^ relative to the Raman crossection of an isolated 
molecule, it varies substantially with nanoparticles shape 
and size. In spheroidal particles, a strong local field en- 
hancement near the tip j2j leads to a lightning rod ef- 
fect recently observed in SERS from nanorods In 
gold nanorods ^2d\ and nanoshells an additional lo- 
cal field enhancement comes from a redshift of the SP en- 
ergy away from the onset of optical transitions between 
electronic rf-band and sp-hand. However, in spherical Au 
particles, the proximity of SP and interband transitions 
energies leads to a damping of SP by interband electron 
excitations and, thus, to a reduction of the SP local field. 



P=aoE 




FIG. 1: Schematic representation of SERS from a molecule 
adsorbed on a metal nanoparticle. 



Such a damping is even stronger in Cu nanoparticles, 
where the SP energy lies above the interband transitions 
onset "2^. In silver particles, however, the interband 
transitions onset lies considerably higher in energy (~ 4.0 
eV) than the SP 3.0 eV for nanoparticles in, e.g., glass 
matrix) and have very litle detrimental effect on SERS. 

Virtually all theoretical studies of the electromagnetic 
mechanism of SERS were performed for relatively large 
nanoparticles with diameters of several dozens nm or 
larger. For such sizes, the SP damping rate, 7, in Ag 
particles comes mainly from the electron-phonon scat- 
tering or electromagnetic retardation effects (for larger 
particles). However, for particle sizes smaller than ~ 10 
nm, the finite-size effects become important. For small 
particles, the width of the resonance peak in absorption 
is determined chiefly by the SP damping due to excita- 
tion of high-energy intraband single-particles transitions 
accompanied by transfer of momentum to nanoparticle 
boundary ■ This effect is usually incorporated via 
size-dependent correction, 7 = 70 -I- 7^ , in the Drude di- 
electric function for sp-electrons. 



s{uj) = 1- ujI/uj{ll> + ij) 



(1) 
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where ujp is the bulk plasmon frequency {cUp ~ 9.0 eV for 
Ag), and 7s is determined by the electron level spacing 
near the Fermi level, 



7s ^ vp/R, 



(2) 



with a numerical coefficient of the order of unity (we 
set h — 1). For small particles, 7s dominates over the 
phonon-induced damping, 70, resulting in a strong de- 
pendence of SERS on nanoparticle size. Indeed, the am- 
plitude of the SP local field at resonance is E oc 7~^, so 
that the electromagnetic enhancement factor |^ depends 
on the radius as A{u}) = \E/Ei\'^ oc (here Ei is the 
incident field). Thus, in small particles, finite size effects 
can reduce SERS by several orders of magnitude. 

In this paper, we demonstrate that for even smaller 
nanometer-sized particles, another quantum-size effect, 
with the opposite trend towards increasing SERS, be- 
comes important. The underlying mechanism is related 
to the difference in the density profiles of sp-band and d- 
band electrons near the nanoparticle boundary. Specifi- 
cally, the localized d-electrons are mainly confined within 
nanoparticle classical volume while the wave-functions 
of delocalized sp-electrons extend outside of it. This 
spillover leads to a larger effective radius for sp-electrons 
[28] and, thus, to the existence of a surface layer with 
diminished d-electron population j2^ 1^ ■ As a re- 
sult, in the surface layer, the screening of sp-electrons 
by d-band electron background is reduced, leading to a 
blueshift of SP absorption peak in small Ag nanopar- 
ticles. More recently, the effect of underscreening of 
Coulomb interactions between sp-electrons has been ob- 
served as an enhancement of electron-electron scatter- 
ing rate measured in ultrafast pump-probe spectroscopy 
[22,1131 3.nd photoemission ^4] experiments. 

Specifically, we address the role underscreening plays 
in the SERS from a molecule located in a close proximity 
to the surface of an Ag nanoparticle. We find a substan- 
tial increase of the SP local field outside of nanoparticle 
and, hence, an additional enhancement of the Raman sig- 
nal. Furthermore, we investigate the size-dependence of 
SERS to find that the enhancement factor deviates con- 
siderably from the behavior when the screening effects 
are included. In particular, the relative enhancement is 
stronger for smaller nanoparticles due to a larger volume 
fraction of the underscreened region. 

The paper is organized as follows. In Section ^ we 
describe the two-region model which we adopt to incor- 
porate the surface layer effect. In Section lllll we calculate 
local fields and Raman enhancement factor. Discussion 
of our numerical results is presented in Section llVI Sec- 
tion concludes the paper. 



II. MODEL 

We consider SERS from a molecule adsorbed on a 
Ag nanoparticle in a medium with dielectric constant 
em- For nanoparticle diameters exceeding w 1.5 nm. 



the bulk electronic structure is essentially preserved while 
the quantum-mechanical corrections due to discreteness 
of the electron energy spectrum can be included via 7s 
in the sp-electron dielectric function Eq. (Q. To incor- 
porate the surface layer effect, we adopt the two- region 
model where sp-band and d-band electrons are confined 
within spherical volumes with different radii Rd and R, 
respectively ^9^, 30, 35]. Note that for d-band elec- 
trons, the semiclassical approach [ssj remains valid even 
for smaller nanometer-sized particles j30j . while for sp- 
electrons, the density ns(r) is a smooth function near the 
boundary while exhibiting Friedel oscillations inside the 
nanoparticles |36|. A fully quantum- mechanical theory 
of SERS in nanometer-sized particles will be published 
elsewhere For nanoparticles under consideration, 

however, the deviations of sp-electron density from clas- 
sical shape do not affect the situation qualitatively p9| 
and ns{r) can be approximated by a step- function with 
a sharp boundary at the effective radius R. 

The frequency-dependent potential is determined from 
the Poisson equation. 



^{uj,r) = Mr) + 



d^r''J^. 



(3) 



where 0o(r) — — eEj • r is the potential of incident light 
with electric field amplitude Ei — EiZ along the ^-axis, 
and 6N{llj, r) is the induced charge density (hereafter we 
suppress the frequency dependence). The latter has four 
contributions, 

(5iV(r) = SNsir) + 6Ndir) + 5N„,ir) + SNo{r), (4) 

originating from the valence sp-electrons, the core d- 
electrons, the dielectric medium, and the molecule, re- 
spectively. The density profile of delocalized sp-electrons 
is not fully inbedded in the background of localized d- 
electrons but extends over it by A = i?— .Rw whfchjjvithin 
our model, is the surface layer thickness ^9, 30, 35]. The 
induced density is expressed via electric polarization vec- 
tor as 

SN{r) - -V ■ P(r) = -V • (P^ + Ps + P,„ + Po), (5) 
with each contribution related back to the potential as 

Pd(r) = -^^^0(i?,-r)V$(r), 



P.(r) 
Pm(r) 



47r 



An 



e{R-r)y^{r), 
r - i?)V$(r), 



Po(r) = -,5(r-ro)aoV$(ro), 



(6) 



where step functions 9{x) enforce the corresponding 
boundary conditions. The molecule is represented by 
a point dipole with polarizability ao located at Fq (we 
chose nanoparticle center as the origin). In the follow- 
ing, averaging over the orientations of molecular dipole 
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is implied so the polarizability tensor ao is assumed to 
be isotropic. 

After substituting above expressions into r.h.s. of Eq. 
P|) and integrating by parts, we obtain a self-consistent 
equation for the potential <&(r), 



e(r)$(r) = 0o(r) - Vo: 



1 



aoVo$(ro) 



+ XseiR-r) + XmO{r' -Rmv'), (7) 

where e(r) — + f-s — ^, Cs, and in the intervals 
r < Rd, Rd < r < R, and r > R, respectively, and 
Xi = (fi ~ l)/47r [i = s, d, m) are the corresponding sus- 
ceptibilities. For simplicity, we assume that the molecule 
is located on the z-axis (along incident field direction). 
After expanding <l>(r) and |r — r'j"-'^ in terms of spherical 
harmonics and retaining only the dipole terms, we obtain 
for the radial component 



the nanoparticle local potential 



1 



Sip{r) 



-l3{r/Rd)MRd 
+Pir/R)MR) 



<Po/e{r) - -E,r/e{r), 
Arf(l — 2A,„) 



1 - 2a^\d\ 
A,„(l - a^Xd) 



1 - 2a^Xd\r. 



(10) 



where a = Rd/R is the "aspect ratio", and the parame- 
ters \d and A„i are given by 



Ari 



ed + 3es - 1 ' 



Am, — 



2e. 



(11) 



The spatial dependence of induced potential, 5<p, is deter- 
mined by I3{x) of Eq. Inside the particle, the poten- 
tial linearly increases for r < Rd , while it exhibits a more 
complicated behavior in the surface layer, Rd < r < R. 
Outside the nanoparticle, 6ip falls off quadratically. 



6(p{r) 



r>R, 



(12) 



e(r)a>(r) = 0o(r) - ^ /3 (^-^) $(i?,) 
47rao ^/ r- \ (9$(ro) 



ra J dro 



(8) 



where 



/3{x) = x-^ e{x - 1) - 2a;6l(l - x). 



The second and third terms in rhs of Eq. © de- 
scribe light scattering from the boundaries at r = Rd and 
r = R, respectively, that separate regions with different 
dielectric functions, while the last term represents the 
potential of the molecular dipole. The boundary values 
of $ are found by setting r — Rd, R,ro in Eq. ((SJ, which 
leads to a closed-form expression for the self-consistent 
potential in the presence of molecule, nanoparticle, and 
dielectric medium. 



III. CALCULATION OF RAMAN SIGNAL 

The dipole moment of a radiating molecular dipole is 
determined by the local electric field, E, at the molecule 
location: p = uqE. The Raman field consists of a direct 
field of this dipole and a secondary field scattered by 
the nanoparticle. In order to extract the Raman signal, 
we present the self-consistent potential in the form $ = 
4> + 4>^, where <j> is the local potential in the absence 
of molecule and cj)^', calculated in the first order in ao, 
determines the Raman signal 0- 

Keeping only zero-order terms in Eq. (|SJ), we find for 



where a{uj) is the particle polarizability 
3 a'^A£i(l — A„i) — A„ 



a(w) = R^ 



1 — 2a'^XdXn 



(13) 



In the absence of the surface layer, Rd — R, we recover 
the usual expression 



id 



1 



ed-l + 2e„ 



(14) 

^=^%=. From Eqs. ^ 

and (|12|l . we then obtain the local electric field at the 
molecule location as 



(9) with the resonance at ujm 



E{ro) 



d(p 
dro 



Ej , , , , a(uj) 

^(l + 2g), g{cu)^^. (15) 



To calculate the Raman signal, we substitute the local 
field into the last term of Eq. |(SJ) and, in the first order 
in ao, obtain the following equation for 0^, 



e(r)0«(r) 



Em, - 1 



ed-l 



f3{r/Rd)cf>"iRd) 



(3{r/R)cf>''iR), 



where the potential 



Ana'r 



^ E{ro) 



2r , 



^eir-ro)~—e{ro 
ro 



(16) 



-(17) 



describes the direct field of a radiating molecular dipole 
(oq is the derivative of molecule polarizabilty with re- 
spect to normal coordinate that determines the Stokes 
shift uJs 113), while the second and third terms describe 
secondary fields due to scattering from d-band and sp- 
band electron distributions boundaries, respectively. The 
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latter can be found by matching (j)^ a.tr = Rd and r = R. 
We now notice that, at these values, (f>l^{r) is linear in r, 
so we can write 0^ = ipo' + 5lpq with 

^^^<P^/eir), S<p^ir)^^Eiro)S^ir), (18) 

where SLp{r) is given by Eq. (|10|1 but with 0^5 instead of 
u!. We then finally obtain for the Raman field (r > tq) 



and, hence, for the enhancement factor 



(19) 



A{lu,lu,) = |1 + 2g{Lo) + 2g{us) + Ag{uj)g{uj,) . (20) 

The above expressions generalize the well-known classical 
result [a S 13 to the case of a small noble- metal particle 
with different profiles of d-band and sp-band densities. 
While SERS retains the usual dependence on nanopar- 
ticle polarizability, the latter is modified in the presence 
of a surface layer [see Eq. ^^]. Note that Eq. ifTHI) 
remains unchanged even for non-classical electron dis- 
tributions provided that electronic wave-functions in a 
nanoparticle do not overlap with molecular orbitals (37i| . 



IV. NUMERICAL RESULTS 

Below we present the results of numerical calculations 
for Ag nanoparticles with diameters ranging from 2 to 6 
nm in a medium with dielectric constant = 2.0. The 
SP resonance is positioned at ujm — 3.0 eV, far away 
from the interband transitions onset in Ag at 4.0 eV, so 
in the frequency range of interest the real part of inter- 
band dielectric function is nearly a constant, ~ 5.2. 
In this size range, the SP damping is dominated by size- 
dependent contribution to 7, with the numerical coeffi- 
cient in 7s [see Eg. (|2Jl] adjusted to fit the experimental 
absorption data |32 |. 

In Fig. 121 we show calculated absorption spectra for 
various surface layer thicknesses A. For finite thick- 
nesses, the SP energy experiences a blueshift whose mag- 
nitude increases with A, in agreement with previous cal- 
culations of absorption in small silver particles 0, [s^l ■ 
This blueshift originates from a reduction of the effective 
(averaged over the volume) interband dielectric function 
in the nanoparticle that determines the SP energy. At 
the same time, the peak amplitude increases with A while 
the resonance width is unchanged. The absolute value of 
polarizability determines, in turn, the magnitude of the 
local field outside the nanoparticle [see Eq. (|15|l ]. 

In Fig. 13 we plot the local electric field at SP reso- 
nance frequency as a function of molecule distance from 
the metal surface, d — r — r^. Outside the nanoparti- 
cle, the local field exhibits the usual decay relative 
to constant incident field background. At the same time, 
the field magnitude is larger for finite A, reaching ~ 20% 



A = 0.0 nm 
A = 0.1 nm 
A = 0.2 nm 
A = 0.3 nm 
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FIG. 2: Absorption spectra of Ag nanoparticle calculated for 
R = 3.0 nm and several values of surface layer thickness A. 



enhancement near the boundary for 3.0 A thick surface 
layer. Enhancement is strongest when the molecule is lo- 
cated in a close proximity to metal surface, near the un- 
derscreened region with small population of d-electrons. 

We now turn to the size-dependence of SERS. In Fig. 
01 we plot the enhancement factor at SP energy as a 
function of particle radius for a small molecule-particle 
separation d = 0.1 nm. We consider here non-resonant 
scattering and assume, as usual, that molecular vibra- 
tional energies lie within SP resonance width (the latter 



A = 0.0 nm 
A = 0.1 nm 
A = 0.2 nm 
A = 0.3 nm 
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FIG. 3: Local field at surface plasmon energy as a function 
of distance from particle boundary calculated for R — 3.0 nm 
and several values of surface layer thickness A. 
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FIG. 4: Size-dependence of Raman enhancement factor cal- 
culated for particle-molecule distance d — 0.1 nm and several 
values of surface layer thickness A. 



being large for nanometer-sized particles). In this case, 
the main contribution to SERS comes from the last term 
of Eq. (|20|) corresponding to the secondary scattered 
field of radiating molecular dipole. As it can be seen in 
Fig. 21 the general tendency is a decrease of the Raman 
signal for small nanoparticles due to a strong SP damp- 
ing by single-particle excitations. For A = 0, the size- 
dependence of enhancement factor is A oc i?'*. However, 
that dependence changes considerably when the effect of 
surface layer is included in the calculations. For finite A, 
the decrease of A is considerably slower, with the signal 
strength in small nanoparticles about 400% larger than 
that for A = 0. It should be emphasized that the role of 
the surface layer becomes more pronounced when parti- 
cle size is decreased. Indeed, the relative enhancement of 
SERS is larger for smaller particles [see Fig. 0] due to a 
larger volume fraction of the underscreened region. 



CONCLUSIONS 



We investigated the role of quantum-size effects in the 
electromagnetic mechanism of surface-enhanced Raman 
scattering in noble-metal nanoparticles. We identified a 
new source of Raman signal enhancement in small par- 
ticles which originates from different density profiles of 
sp-band and d-band electrons near the boundary. The 
existence of an underscreened surface layer with low pop- 
ulation of d-electrons gives rise to a stronger, compared 
to classical calculations, local field of surface plasmon col- 
lective excitation that determines the magnitude of Ra- 
man signal from a molecule near the surface. Although 
the dominant finite-size effect is still a reduction of SERS 
due to surface plasmon damping, the decrease of the Ra- 
man signal is considerably slower when the surface layer 
effect is taken into account. 

We calculated the size-dependence of SERS in the 
framework of two region model with different (but classi- 
cal) distributions of sp-band and d-band electron densi- 
ties and semiclassical treatment of electron surface scat- 
tering, and found that the additional enhancement be- 
comes stronger as the particle radius decreases. Although 
for particles size smaller than 2-3 nm, the semiclassical 
model is no longer valid, the physical mechanism of the 
enhancement persists for even smaller sizes. In fact, the 
semiclassical model understimates the fraction of under- 
screened region as the particle size decreases. On the 
other hand, in the absence of a sharp boundary the lo- 
cal fields are weaker. The outcome of this competition 
depends on the precise shape of the electron density as 
well as on the surrounding dielectric. A fully quantum- 
mechanical calculations of SERS will be presented in a 
subsequent publication Is^ . 
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